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Abstract 

In this work, we study the Brinkman-Forchheimer equations driven under shp 
boundary conditions of friction type. We prove the existence and unique- 
ness of weak solutions by means of regularization combined with the Faedo- 
Galerkin approach. Next we discuss the continuity of the solution with respect 
to Brinkman's and Forchheimer's coefficients. Finally, we show that the weak 
solution of the corresponding stationary problem is stable. 

Keywords: Brinkman-Forchheimer equations. Slip boundary conditions. Weak 
solutions. Continuous dependence, Stability. 

AMS Subject classification: 35J85, 35Q30, 76D03, 76D07 

1 Introduction 

We consider the Brinkman-Forchheimer equations for unsteady fiows of incom- 
pressible fiuids, i.e. 

du 

— -lyAu + au + blul^u + Vp = / in Q = O x (0, T) , (1.1) 
at 

divu = in g, (1.2) 

where fl is the fiow region, a bounded domain in . The motion of our in- 
compressible fiuid is described by the velocity u{x,t) and pressure p{x,t). In 
(jl.ll) and (|1.2[) . / is the external body force per unit volume depending on x 
and t, while the positive parameters v, a, b are respectively the Brinkman coeffi- 
cient, the Darcy coefficient and Forchheimer coefficient, and a G [1, 2] is a given 
number. Equations (jl.ip and (|1.2p are supplemented by boundary and initial 
conditions. As far as the initial condition goes, we assume that 

M(-,0)=tto on , (1.3) 

where no is a given function, that will be made precise later, and fl is the 
closure of ft. Next in order to describe the motion of the fluid at the boundary, 
we assume that the boundary of fl, say. Oil, is made of two components S (say 
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the outer wall) and F (the inner wall), and it is required that d^l = S UT, with 
5 n r = 0. We assume the homogeneous Dirichlet condition on F, that is 

M = onrx(0,r). (1.4) 

We have chosen to work with homogeneous condition on the velocity in order 
to avoid the technical arguments linked to the Hopf lemma (see [15] , Chapter 
4, Lemma 2.3) . On S, we first assume the impermeability condition 

UN = u-n^Q on S X {0,T) , (1.5) 

where n is the outward unit normal on the boundary dfl, and un is the normal 
component of the velocity, while Uj- ~ u ^ unti is its tangential component. 
In addition to (|1.5p we also impose on 5", a threshold slip condition [211 [7], 
which is the main ingredient of this work. The threshold slip condition can be 
formulated with the knowledge of a positive function g : S — > (0, oo) which is 
called the barrier of threshold function and the use of sub-differential to link 
quantities of interest. It is written as 

- {crn).r e gd\u.r\ onS'x(0,T), (1.6) 

where (crn),- is the tangential component of the Cauchy tensor <t given by 
a- = -pI + 2vD{u) with D{u) = ^[Vu+ (Vw)^], and 9| • | is the sub-differential 
of the real valued function | • |, with \w\'^ — w ■ w . We recall that if X is a 
Hilbert space with xq € X , then 

y G 9'I'(xo) if and only if ^"(2;) - *(a;o) > 2/ • (.T - a:o) \/x e X . (1.7) 

Without using the sub-differential, the threshold condition (|1.6p can be written 
as [5] 

|(cr'^)T| < g, 
Kern),-! < g ^ u-r = 0, 

Kern),- 1 = g^u-r ^0 , 

One observes that different boundary conditions describe different physical phe- 
nomena. In [31] , the equations of Brinkman corresponding to (|1.1|) with 6 = 
have been derived using mixtures theory, in fact a class of approximate mod- 
els for flows of incompressible fluids passing porous solids have been described. 
Forchheimer [6] studied flow experiments in sandpacks and came to the conclu- 
sion that for small Reynolds numbers (Re « 0.2), the diffusion law of Darcy 
is not significant. Furthermore, he found the relationship between the pressure 
gradient and the velocity obtained using the law of Darcy to be nonlinear. In 
fact for a wide range of physical experiments, he found that the nonlinear term 
should be quadratic. Inertial effects in the porous medium at relatively small 
Reynolds numbers are the cause of the nonlinear excess pressure drop observed 
by Forchheimer and others. The slip boundary conditions of friction type ()1.6p 



on S'x(0,T). (1.8) 
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can be justified by the fact that frictional effects of the fluid at the pores of 
the sohd can be very important. In fact on the role of the boundary condi- 
tions for such problems, Brinkman [3] mentioned that "The flow through this 
porous media is described by a modification of Darcy's equation. Such modi- 
fication was necessary to obtain consistent boundary conditions" . While there 
continues to be some debate over the functionality of the Brinkman- Forchheimer 
model [28j , nonlinearity has been verified experimentally (22j , and some theoret- 
ical results have been obtained in [H [521 HI E 130] • The Brinkman- Forchheimer 
equation continues to be used for predicting the velocity of the flow in the 
porous region, while the energy equation for the porous region accounts for the 
effect of thermal dispersion [53]. Since we are well aware that for such flow, 
there are important features at the boundary, it is therefore important to model 
Brinkman-Forchheimer flow accurately taking into account the motion at the 
boundary. This is the driving force behind our work. 

Even though flows under boundary conditions of friction type have been consid- 
ered in various publications ([Ml[Zl[lSl[Il[ni[Il[Hl[ni[Tn] among others), and 
Brinkman-Forchheimer equations (|l.ip , (|1.2p with non slip boundary conditions 
has been examined in [1 [29l [2 [T] [30] , the combination of (frTj) . (fT^ and pTH]) 
has not been presented in the literature, and it is the object of this work. The 
novelty of the problem, from the mathematical point of view, derives from the 
following features; the highly coupled and nonlinear nature of the problem, the 
incompressibility constraint and related pressure, and the leak boundary condi- 
tions ([13]) and ([L6)) . 

Not surprisingly, flows problems involving boundary conditions of friction type 
offer significant theoretical and computational challenges. With regard to theo- 
retical studies, the work by Hiroshi Fujita and co-authors [7] [T3 l [TT ] [T2 l IS] [9] [TO] . 
represent some early, contributions. These authors established existence, and 
uniqueness of solutions, for the equations corresponding to Stokes equations 
by means of semi-group approach, regularity of solutions arc also examined. 
An interesting and related work is that by Christiaan Leroux and co-author 
[Ml [55] on Stokes and Navier Stokes equations under more general "friction 
type boundary conditions". As far as computational studies for flows driven 
by "friction type boundary conditions" are concerned, it should be mentioned 
that even though the literature is very rich in problems formulated in terms of 
variational inequalities [171 [El UHl US : not much have been done for the specific 
case involving mixed coupled problems [HI [26] [HI [33] [21] , and we would like to 
explore that research direction. 

Problem (|l.ip - (jl.6p is a coupled nonlinear system of equations with a non- 
differentiable expression (at zero) on S due to the sub-differential term i9|mt|- 
We propose to solve the resulting system of partial differential equations us- 
ing the regularization approach [S] I27j , which consists of replacing the original 
problem by a sequence of "better behaved" approximate problems indexed by 
a small positive parameter e. We then solve the regularized problems by the 
Faedo-Galerkin method, and finally, the solution of the original problem is ob- 
tained by passage to the limit as e goes to zero. The difficulty in the algorithm 
described is to obtain the pressure. Indeed, as the problem in its weak form 
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is formulated as a variational inequality with only one unknown, the pressure 
will not be obtained in the usual way (for the classical Navier-Stokes equations 
see e.g [32], [Theorem 2.5-1, page 54]). But, instead we first construct a regu- 
larized pressure by using the classical approach and then pass to the limit as e 
goes to zero, after showing that the regularized pressures are bounded in some 
appropriate function space. After constructing weak solutions of the problem, 
we analyze some qualitative properties of the solution, namely; the continuous 
dependence of the solution with respect to the Brinkman and Forchheimer co- 
efficients, and the stability of the stationary solution. The results presented, 
extend in some sense those obtained in [521 ISO] to a family of variational in- 
equalities with non-differentiable functionals. 

The remaining part if this work is organized as follows. In section 2, we docu- 
ment the variational formulation associated to the problem and prove its well- 
posedness. Section 3 is devoted to the stability of the solutions with respect to 
some data of the problem. The stability of the stationary solutions is analyzed 
in Section 4. 

2 Analysis of the problem: Solvability 

We introduce some preliminaries and notation for the mathematical setting of 
the problem. We write down a variational formulation of problem (jl.ip - (jl.6p . 
Next wc derive some a priori estimates of its solution and obtain the existence 
of solutions by means of Faedo-Galerkin. 

2.1 Preliminaries/Notation 

In what follows, for 1 < p < oo. LP(f7), and LP{d^l) are the usual Lebesgue 
spaces, with norms denoted by || • jj^p and || • ||LP(asi) respectively, (of course when 
p = 2, we will denoted the norm in L^{^1) by || ■ ||). We shall use the following 
notation; for the sake of simplicity, one defines them in three dimensions. Let 
k = (fci, k2, ks) denote a triple of non-negative intergers, set |fc| ~ ki + k2 + k^ 
and define he partial derivative d'' by 



Then, for non-negative integer m, we recall the classical Sobolev space 



d\k\y 



H"\n) = {v e L^{n) ; d''v e L^{Vl) V |fc| < m} 



equipped with the seminorm 




1 1/2 
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and norm 



1 1/2 



0<fe<Tn 



For p = 1,2,3, the inner products in the spaces L'^{d^Y and 

H^{^Y are denoted by (•, •), (■, and (•, •)!, respectively. The product spaces 
are denoted by bold letters: VL^U) = H\Vlf , I? in) = L^iUf, L"+2(rj) = 
etc. 



Here, and in what follows, the boundary values are to be understood in the 
sense of traces. We omit the trace operators where the meaning is direct; oth- 
erwise we denote the traces by z;|r, v\s, etc. Also, all the derivatives should be 
understood in the sense of distribution. 

Wc also recall from [15] (Chap. I, Thm 1.1) for instance the following Poincare- 
Friedrichs inequality: 



for all t;eHi(r!)nK|s = 0, t;|r = 0}, ||-u|| < C||Vt;|| , (2.1) 



which yields the equivalence of the norms || • |ji and | ■ |i on H^(f2) n {v„|s = 
, t;|r -0}. 

For any separable Banach space E equipped with the norm || ■ H^, we denote 
by C"(0,r;i?) the space of continuous functions from [0,T] with values in E 
and by D' {0,T; E) the space of distributions with values in E. LP{0,T; E) is 
a Banach space consisting of (classes of) functions t i — > f{t) measurable from 
[0, T] I — > E (for the measure dt) such that 



In what follows, (f){t) stands for the function x ^ il i-^ 4'{x,t). 

We assume that the data {f,g) belong to L'^{0,T;L'^{n)) x L°°{SY, and that 

the datum Mq belongs to H^(fi) n L"+^(51), and satisfies the incompressibility 

condition 



This last condition is not necessary for all the results that follow but, since it is 
not restrictive, we shall assume it from now on. 

2.2 Variational formulation 

In order to write a variational form associated with (|l.ip - (|1.6p . we retain ()1.3p 
and we weaken the equations (|1.1|) . (|1.2|) and constraints (|1.4|) . (jl.Sp using the 
Green's formula, while (|1.6p is re- interpreted with the help of (|1.7p . It follows 
from the nonlinear term in (|l.ip that u{t) and the test function v should be- 
long to L"+^(r2). Then u'{t) and \u{t)\'^u{t) must belong to the conjugate of 
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L~(0,T;_E) = eSSo<t<TSUp ||/(t)||£; < CO . 



div Mq = in 51 . 



(2.2) 
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L"+^(ri). which is L"+i (ft) . Wc then introduce the following spaces 

JV = Hi(r!)n{-y|r =0 , i'„|s = 0} , 
M = Ll{n) = {q e L^{n) , {q,l) ^0} . 

We then adopt the following definition of weak solutions of (|l.ip - (|1.6p 

Definition 2.1 Given {f,g) m L^{0,T;L^{fl)) x L°"{S)^, and Uq € U^ifl) n 
L""'"^(n), satisfying \2.2jl . We say that {u,p) is a weak solution of li. J|) - |77^j 
if and only if; 

u€L°° (0, T; {N n L"+2 {n)),pcL^{Q, T; M) , and u' G (O, T; {N n L"+2 (f]))') , 
and 

for almost all t and all q e L'^{^),v E M n L"+2(ri) 

iu'{t), V ~ u{t)) + 7(Vit(t), V{v - u{t)) ) + a iu{t),v - tt(t)) 

+ 6 V - u{t)) - (div(t> - u{t)),p{t)) + 

+Jiv) - J{u{t)) > {f{t),v - u{t)) , (2.3) 

(divM(<),g) = (2.4) 

u{0) = uo , (2.5) 

where, J{v) {g{x), \vt{x)\)l2(s) ■ 

Following [5], it can be shown that any solution of (jl.l|) - (|1.6p is a solution 
of (|2.3p - (|2.5p in the sense of distributions. The converse property holds for 
any solution of the problem (|l.ip - (|1.6p that enjoys the regularity mentioned in 
Definition 12. 1[ in a sense to be made precise later on. The kernel of the bilinear 
and continuous form L'^{D,) x TV n L"+2(fJ) 3 {q,v) i — > {q,divv) e M is 
Y = {v € Afn L"+2(f]) , div t> = in n} . With the space V in mind, it is 
then easy to see that the function u{t) given in p.Sp - p.Sp is a solution of the 
simpler variational problem: Find it e L°°{0, T; Y) , u' e L^(0, T; V') satisfying 
(|2.5p such that for almost all t and all ?; G V 

{u'{t), V - u{t)) + j{Vu{t), \/{v - u{t)) ) + a {u{t),v - u{t)) (2.6) 
+b{\uit)\"u{t),v - u{t)) + Jiv) - J{u{t)) > if{t),v - u{t)) . 

Next, we establish the solvability of the variational problem (|2.6p by means of 
regularization combined with Galerkin's method. We then construct a pressure 
p in L^(0, T; ig(il)) such that the couple {u,p) enjoys the regularity announced 
in definition 12. 1[ and satisfies (|2.3p - (|2.5p . 

2.3 Existence of a solution 

In this paragraph wc discuss the solvability of ()2.6p by regularization, and pas- 
sage to the limit. Thus it is obtained in several steps, that we describe below. 

Step 1: Regularized problem. 
We first recall that one of the difficulties of solving (|2.6p is the fact that the 
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functional v £ V i — > J{v) = {g{x), \vt{x)\)s is not differentiable at zero. So, 
to bypass that hurdle we introduce the regularized functional defined by 

V <eY ^ Je{v) = {g{x), ^/\vt{x)\'^ + e'^)s , < e << 1 . 

Clearly is convex and Gateaux differentiable with Gateaux derivative Kg 
defined on V and given by 

{KJu),v) ~ / q — =ds. 

Wc briefly observe that K,, is monotone. Indeed since is convex, for m, v 
elements of V and < t < 1, Je(tM + (1 - t)v) < tJ^{u) + (1 - t)J^(v), which 
can be re- written as 

Je{v + t{u-v)) -J,{V) , J 

Then by taking the limit on both sides when t goes to zero yields 

(A'e(t)), U-v)< J^{u) - J^{v) . 

Interchanging the role of v and u, one gets instead 

{K,{u),v-u)<J,{v)-J,{u). 

Putting together the later and former inequality, one has the desired result 

{K,{u) ^ Kg{v),u~v) >0 . (2.7) 

The regularized form of (|2.6p can be written as follows: Find Me G L°°{0,T;Y) 
satisfying (|2.5p with u'^ e i^(0,r; V') such that for almost all t and aU t) £ V 

+bi\u,it)\^u,it),v - u,it)) + Mv) - Je{u,{t)) > (fit), V - u,it)l2.8) 

Since Jg is differentiable, adopting the classical arumcnts in [5], one can state 
that ([278]) is equivalent to: Find G L°°(0,T;V) satisfying (|23)) with u'^ G 
2.2(0, T;V') such that for almost all t 

{u'^{t), v) + u{Vu,{t),\/v) + a {u,{t),v) + b {\u,{t)\''u,{t),v) 

+ {K,{u,it)),v) ^{fit),v) forall-uGV. (2.9) 

Before proving the existence of a solution ^^(t) of (|2.9p . wc first show how the 
pressure is constructed, knowing the velocity. For that purpose, we begin by 
integrating (HH) on [0,t], apply ([23]), and for v G A/" n L"+2(rj); we introduce 
the functional 

Hiv)it) = f [(/(s), v) - u{Vu,{s),Vv) - a (m,(s), v) - b i\u,is)ru,is),v) 
Jo '- 

t 

{Kg{ug{s)),v) ds - {u.{t),v) + {uo,v) , for aU 0<t<T. 
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One sees that 'H is linear and continuous on A/" n L"+^(r2), and according to 
(|2J)) and ((23)) . it vanishes on V. Now following [32], [Theorem 2.5-1, page 54], 
for each t S [0,T], there exists a unique function Pe(t) G io(^) ^ positive 
constant C such that: for all v e A/" n L"+2(rj), 

HW(0 = (divi;,p,(0) , (2.10) 
C|b~ (Oil < sup (-^^-;-f . (2.11) 

Finally, we take the time derivative on both sides of (|2.10p : and we let 

pS) = j^Pe{t) , (2.12) 

in the resulting equation. Thus we have obtained the following variational 
problem: Find € £^(0, T; A/" n L"+2(17)), p. e L'^{0,T; Ll{n)) with u'^ G 
L2(0,T;(A/'n L"+2(f2))') such that for almost all t and ah q e L^{n) , v G 
A/'nL"+2(17) 

(u^(<), v) + iy{Vu,{t),\7v) + a {u,{t),v) + b {\u,{t)\''u,{t),v) 
-idWv,p,{t)) + {K,{u,{t)),v) = {f{t),v), (2.13) 
(divMe(t),g) = 

Me(0) = Mo . 

It is clear that the variational problems p.9p and p.lSp are equivalent, with 
the regularized pressure described by p.lOp , (|2.1ip and p.l2p . 

Step 2: Faedo-Galerkin approximation. 
Wc let 

H = {veL'^in), divt; = 0, w„|ao =0}nL"+2(17) . 

One readily observes that V is compactly embedded in H . For the slip bound- 
ary condition, we introduce the Stokes operator defined on a subspace of V 
constructed in [20] as follows; for every / g H, there exists a unique v such 
that 

(Vi;,V(^) = (/,0),V0eV. (2.14) 

Moreover, for every d G V, there exists a unique / G HI such that (|2.14l) holds. 
Then p.l4p defines a one-to-one mapping between / G H and v G D{A), 
where D{A) is a subspace of V. Hence, Av = f defines the Stokes operator 
A : D{A) — ^ H. Its inverse A^^ is compact and self-adjoint as a mapping from 
IH to H and possesses an orthogonal sequence of eigenfunctions ipf. which are 
complete in H and V; 

ArP^ = AfcVfc. (2.15) 
Let Vm be the subspace of V spanned by -01, • . • , ipm, that is 
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We consider the following ordinary differential equation: Find Ug^„i{t) £ V,„ such 
that for all v G 

(u^ .^{t) , v) + v{Vus^„i{t) , \/v) + a{u^,,n{t) , v) 
+b{\u,,n{tTu,,n{t) , v) + {K,{u,,n{t)) , v) = {f{t),v) , (2.16) 

Us,,niO) Me(0) ^ Uq £ V™ . 

As far as the existence of u^^m{t) defined by (|2.16p is concerned, we note that 
the mapping 

K, : w I — > (/, v) — v{'\Iw , \7v) — a{w , v) — b {Iwl^w , v) — {K^{w) , v) , 

is locally Lipschitz thanks to the nature of the operators involved. It then follows 
from the theory of ordinary differential equations that (|2.16p has a solution M^.m 
defined on [0,te,,n], t^.m > 0. Hereafter, C denotes a constant independent of 
m, and depending only on the data such as fi, and whose value may be different 
in each inequality. Next, we derive some a priori estimates and deduce that 
te^m does not depend on s or m. Concerning the later property, is should be 
mentioned from [SJ I27j , that it suffice to derive a priori estimates of the solution 
with the right hand side independent of m and e. 

Step 3: Some a priori estimates. 
First we let v — u^,rn{t) in (|2.16|) . After using Young's inequality, one obtains 



d_ 

dt 



iMe.mWf + 2i^||VM,,™(t)f + a||tt,,,„(t)f + 2b\\u,^„,{t)\\ 



\\f(t)\\^ 

+2(i^,(M,,,„(t)) , M,.,„(i)) < ^ (2.17) 

a 

which by integrating over [0,r»] for r» < te,m, and using (|2.7p . yields 

sup ||Me.r„(t)||' + 2i. / 1 1 Vm,,„ (f) | j^dt + / | (t) 1 1 ^di 

0<i<T9 Jo JQ 

+2b [ \\u,,,nmitl2dt<- [ \\f{t)\\^dt + \\uo\\^ <oo (2.18) 
Jo aj„ 

since by assumption / e L'^iQ) ■ Now let v = u^^^{t) in ^IJ^. For < t < T^*, 
Young's inequality yields 



-|[2/Ve(«,.™(i))]<||/(t)|P 



which leads to 



26 



\\u^,„M\'dt + y\\Vu,^^{t)\\' + a||M,,™(t)||^ + -^\\u,,^{t)\\ltl. 

a + 2 

2K,{u,.ra{t)) < \\f\\l2^Q^j.t^2^+m, (2-19) 
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where 

$(0) = zy||VMo|P + a||Mo|P + —^\\uo\\ttl2 + 2 / g^/\uo\■^ + l ds . 

It is manifest that the right hand sides of the a priori estimates obtained in 
(|2.18p and (|2.19p are independent of m and e. We then conclude that te^m is 
independent of e and m foUowing the arguments discussed in length by [SJ [37] . 

Step 4: Passage to the limit. 
Wc need to pass to the limit when m approaches infinity and e approaches zero. 
We start by fixing e and study the sequence m i — > M£,m- 
Based on (|2.18p and (|2.19p . it is clear that when m — > oo, 

u^^rri remains bounded in L°° {0,T;M), 

|Me.„i|"Me,m remains bounded in (0, T; (n)) , (2.20) 
^ remains bounded in L^{0,T;L^{n)) . 

From a consequence of the result of Dunford- Pettis [33] , it is possible to extract 
from a subsequence, denoted again by (we^m)™ such that 

Ue^m — > Ue weak star in L°°{0, T; H) (2.21) 
Me,™ — > Me weak star in L°°(0,T;V,„) (2.22) 

— > Xe weak star in ( 0, T; L°+i (fi) 1 (2.23) 

„ — > weak in L^{Q,T;m). (2.24) 

The convergence results p.2ip . and p.24p imply in particular that 

Ue.m remams m a bounded set of [Q] . (2.25) 

But from Rellich-KondrachofF, the embedding H^{Q) \ — > L^{Q) is compact. So 
one can extract from (Me,m) a subsequence, denoted again by (ite.m) such that 

Ue.m — > strong in L^(0,T;H) and a.c. in Q . (2.26) 

Next, it follows from and and Lemma 1.3 in gT] (page 12) that 

It remains to be shown that 

Keiue,rn) — > K^{u^) wcak Star in L°°(0, T, V^J . (2.27) 

Firstly from ([^T^ 

Keiue^ra) — > Pe weak star in L°°(0,r,V;j . (2.28) 

Passing to the limit in (|2.16p . one obtains 

{u'^{t),v) + v{Vu,{t) , Vv)+a{u,{t) , v) + b{\u,{t)\°'u^{t) , v) 

+ {l3,,v)^{f{t),v),yveY„. ■ (2.29) 



10 



For any w G L^{0,T;Y„i), since i^e(-) is monotone fsee l2.7p . 
but from ((2l6ll 

{Ke{Ue,m{t)) , U^^rn{t)) = {f {t) , Ui^^rnit)) - {u^„^{t) , 

-:/(VMe_,„(t) , VMe,„i(t)) - a (t) , Me,„i(t)) 



The former and latter equations give 
1 d 

> {Ke{Ue,m{t)) , w) + {R'eiw) , Ue,m{t) - w) 



So, integrating with respect to t on [0, T], yields 



'^(/(0,Me.™(0)rft-^ll«e,rn(r)f + i||M,,„,(0)f 



[l^\\Vu,,,nit)f + a\\UeMt)f +HUe,r4t)rLtl2] dt (2.30) 
> / [{K,{u,^m{t)) , wit)) + {K,{w{t)) , Me,m(t) - 11^(0)] ■ 

Jo 

Next, wc take v ~ Ue,m{t) in p.29p . and combined the resulting equation with 
(|2.30|) . which yields (after taking the limit as m approaches to infinity) 

T 

(/3e - KAw{t)) , u,{t) - w{t))dt > . (2.31) 

At this juncture, we let u^it) - w{t) = ±q with q e L'^{0,T-Ym) ■ Thus (|23T|) 
leads to 



/ {13, - K,{w{t)) , q)dt = , 
Jo 



from which we deduce the desired convergence result (|2.27p . We have established 
that as m goes to infinity, the sequence (we,m(i))m converges to Ug^t) in some 
sense with Ug{t), the solution of 

«(t),7;) + i.(VM,(0 , Vv) + aiu,it) , v) + b i\u,{t)\''u,it) , i;) 

+ (if,(M,(t)),-(;) = (/(i),t;) , for all £ V™ . ^ 

Since U„iV„i is dense in V, we can conclude that (|2.32p holds true for v 
in V. Therefore, we have established that there exists a function uniformly 
bounded with respect to e in L°°(0, T, ]HInVnL"+2(fl)) such that u', is uniformly 
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bounded with respect to e in L^(0,T, H) and satisfies (|2.32p . 

Our final task in the paragraph is to consider the limit as e goes to zero. 

First, we take the limit on both sides of (|2.18p and (|2.19|) . one has 



sup + 2i^ / \\\7ue{t)\\'^dt + a [ \\ue{t)\\'^dt 

0<t<T Jo Jo 

+26 r \\ue{t)r+Ut<- r \\f{t)\\^dt+\\uo\\', (2.33) 
Jo a Jo 



and 



T 



K(t)|pdi + H|V«,(t)|P+a||«,(i)ir + ^||«e(t)ll^^^. 

a + 2 



< I ||/(t)||2dt + $(0) . (2.34) 
Jo 

Thus we can extract from (Me)e a subsequence still denoted by {ue)e such that 

Me — >u weak star in i°°(0, T, H) (2.35) 
Ue — >u weak star in L°°(0, T, V) (2.36) 

\ue\°'Ue — > X weak star in (^0,T,L^(rj)^ (2.37) 

— > u' weak in L^{Q,T,m) . (2.38) 

Arguing as before we can prove that 

Ue — > u strong in L(0,r;H) and a.e. in Q, (2.39) 

\ue\°'Ue — > weak in (^0,r;L^(f7)) . (2.40) 

Let V e £2(0 r from (|232|) . it follows that 

{u'.St).V-Ue{t)) + y{yUe{t) , V [v - Ueii))) + a {Ue{t) , V - Ue{t)) 
+b{\Ue{t)\"Ue{t) , V - Ue{t)) + Je{v) - Je{Ue{t)) 

= (/(t), - Ue(t)) + J,(t;) - Je{Ue{t)) - {Ke{Ue{t)), V - Ue{t)) . (2.41) 

Integrating (|2.4ip with respect to t along [0,T] and taking into account the fact 
that Je{v) — Je{ue{t)) — {Ke{ue{t)),v — Ue{t)) > 0, onc obtains 

T 

{{u'^{t),v) + !^{VUe{t),\/v) + a{Ue{t),v) + b{\Ue{t)\"Ue{t),v)) dt 

+ [ {Je{v)-{f{t),V^Uemdt 

J, (2.42) 

> IwMTW - ^ll^^olP +^ (a||«s(i)lP + \ue{tr+'dx^ dt 

+ f Je(,Ue{t))dt . 
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Since — u weak star in L°°(0,T, V), and Jg is a convex and continuous 
functional on V, one has 

lim inf / J.[u,{t))dt> I J[u{t))dt . (2.43) 
^-^0 Jo Jo 

By using ([2^ . we infer from (|2^ that 

{{u\t), v) + iy{Wu{t), Vv) + a{u{t), v) + b{\u{t)\''u{t),v) 

+J{v)-{f{t),v-u{t)))dt 
> i||M(r)|p-i||Moir+^ (^a\\u{t)\\^ + b \u{tT+^dx^dt + J{u{t))dt 

[(u'it), uit)) + a iu{t),u{t)) + b {\u{t)\"uit), uit)) + Jiu{t))] dt 

which by arguing as in [5], pages 56-57, yields 

{u'{t),v - u{t)) + v{Vu{t),V{v - u{t))) + a{u{t), v - u{t)) 
+b{\u{t)\'^u{t),v-u{t)) + J{v)- J{u{t))>{f{t),v-u{t)) for aU v eY . 

We then conclude that 



Theorem 2.1 The variational problem S2. 9\) admits at least a weak solution, 
which moreover satisfies; 



sup ||VM(t)||<C , / \\u'{t)\\^dt<C 

0<t<T Jo 



(2.44) 



where C is a positive constant depending on the data. 



Having obtained the velocity, we shall indicate how the pressure is constructed. 
First, we recall that from (|2.13P i . 

{divv, ps{t)) = {u'^{t),v) + iy{\7u.{t),\7v) +a{Us{t),v) 
+b{\u,{t)ru,{t),v) + {K,{u,{t)),v) - [f{t),v) , 

but since Pe{t) G L'^{Vl), following [TS], one can find a positive constant C such 
that 

n\\ ^ (divt>,pe(t)) 

C\\pe{t)\\ < sup — . 

veAf \\v\\i 

Now, combining the former and latter equations and the continuity of operators 
involved, one obtains 

C\\p,{t)\\ < ||<(t)||+i.|IVMai)ll+«lhs(i)ll+&ll«e(t)|l2tV. 

+ \\K,M))\W, + \\fit)\\ 
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Equivalently, 



la+l 



CWpeim < ||w:.(t)|| +i^||VM,(t)|| +a||MeWII +C(6,l},a)|iw,(i)llla 
+ Cmgh^.^s)\\Mt)h + \\fit)l 
< \\u'^it)\\ + i^\\\7u,{t)\\ + a\\u,{t)\\ + C{b, n,a)\\Vu,{t)r+^ 

+ Cm\9\\L^iS)\\Ue{t)h + \\fm 

which by Young's inequaUty and integrating the resulting inequahty over [0, T], 
yields (after utilization of ((OH)) and ((011) ) 



\\MtWdt<C \\u'Mf + C \\Vu,it)fdt + C \\u,{t)f 
Jo Jo Jo 

+C r \\Vu,it)r'^+'dt + C\\g\\l^^s) f \\^e{t)\\l 
Jo Jo 

+c r \\f{t)fdt<^ , 

Jo 

(2.45) 

C being a positive constant depending on the parameters and the domain of 
the problem. Then we can select from Pe{t) a sequence, again denoted by Pe(i), 
such that 

p, — >p weakly in L^{0,T; Ll{n)) . (2.46) 

Next, one observes that p.l3p can be re- written as 

{u'^{t),v - u,{t)) + i/(VM,(t), V(-i; - u,{t))) + a {u,{t),v - u,{t)) 
+b{\u,{t)ru,{t),v - U.S)) - (diviv - u,{t)) , p,{t)) 
+Jeiv) - J^iueit)) - {f{t),v -u,{t)) >0 for all t> e A/'nL"+2(17), 

(divMe(i),g) = ,for ah qeL^{n), 

which by integration over the time interval [0,r] and passage to the limit (as 
e ) yields, (after utilization of the identity (div u^{t),q) — for all q G 

T 

[{u'{t),v - u{t)) + v{Vu{t),S/{v - u{t))) + a {u{t),v - u{t))] dt 

+ [ [- {dW{v - u{t)) , p{t)) + J{v) - J{u{t)) - (fit), V - u{t))] dt 
Jo 

+ b{\u{t)\"u{t),v-u{t))dt>0, 
Jo 

for all V eJVn L"+2(rj). Also, (divM(i), q) = for all q G L^^n) . 



14 



Finally, arguing as in [S], pages 56-57, one obtains 

{u'{t), V - u{t)) + v{Vu{t), V{v - u{t))) + a {u{t), v - u{t)) 
+b {\u{t)\°'u{t),v -u{t)) - {div{v -u{t)) , p{t)) + J{v)-J{u{t)) (2.47) 

>{f{t),v-u(t)) 

for all t> e A/" n L"+2(17). Moreover, (div u{t),q) = for all q e L^{n). 

3 Continuous dependence on the data 

In this section, our focus is to establish some qualitative properties of the weak 
solutions in Theorem 12.11 In particular, we show that the solutions depend 
continuously on initial velocity, external force as well as the Forchheimer's and 
Brinkman's coefficients. We recall that such results in the literature are some- 
times referred to as structural stability. 
We first claim that 

Theorem 3.1 Let Ui he the solution of V2. 5]) with respect to UiQ,f^, i = 1,2. 
Then there exists a positive constant C, depending on a, v and J7 such that 

- M2(<)IP < e-'^*|lwi(0) - uM\? + r e'='^"*+'^|I/2(s) - SM\?ds. 

Jo 

(3.1) 

This theorem implies in particular the following uniqueness result. 
Corollary 3.1 The problem h2. 5\] has one and only one solution. 
Proof theorem 13.11 The functions Ui and U2 satisfy respectively: 

{dtUi,v — Ui) + v{yui,V{v — Ml)) + a {ui,v — Ui) + b {\ui\"ui,v — Ui) 

+J{v) ~ J{ui) > (/i,i; - Ml) for aU v eY . 

(3.2) 

and 

{dtU2, V ~U2) + v{yU2, V{V - U2)) + a {U2,V - U2) + b {\U2\°'U2, V - U2) 

+J{v) - J{u2) > {f2,v- U2) for all v eY . 

(3.3) 

Setting v = tt2 in p.2p and t; = iti in (|3.3p and adding the resulting inequalities, 
it follows that 

^±\\wit)\\^ + i.\\S7w\\^ + a\\wit)\\^ 

+ bi\u2rU2 - \ui\'-Ui,w{t)) < if 2 - fi,w{t)) , 
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where w{t) = it2(t) — tti(i) and Wq = tt2o — Mio- Since T(C) = |C|"C is monotone 
then 

i\u2\''U2-\ui\''Ui,w{t))>0 . 

Therefore 

j^\\w{t)\\' + Ciiy,a,n)\\w{t)\\' <C{a,n)\\f2- , (3.4) 

where Poincare's inequahty has been used. We readily deduce the desired result 
from (|3.4p using Gronwall's lemma . □ 

In line of theorem 13. li one can state the following result. 

Theorem 3.2 The weak solutions of problem 112. 5]) constructed in theorem \2.1\ 
depends continuously with respect to the norm on: 

(a) the Forchheimer coefficient h, and 

(b ) the Brinkman coefficient v . 

The proof follows mutatis mutandis the proof of theorem 13. II . 

4 Stability of stationary solutions 

Hereafter, we study the stabihty of stationary solutions to (|2.5|) . 

We assume that the apply force / is independent of time, and we consider the 

following stationary problem 

-v^u + ait + 5jtt|"M — Vp = /, in VL , 
div M = 0, in 17 , , ^ 

It = on 1 , 

tt • n = , and — ctt E gdlu-r] on S . 

Here, we always assume that a G [1, 2], "/,a,b > 0. 

It is clear that the velocity satisfies the simpler variational problem 

{Find It e V such that for all t> G V , 
i^(Vtt, V(f — u)) + a {u, V — u) + b{\u\"u, v — u) + J{v) — J{u) > (/, v — u) 

(4.2) 

It can be shown as in [5] that there exists a unique m G V such that (|4.2p holds 
true, and one has the following 



Theorem 4.1 The weak solution u of S2. 5\) constructed in theorem \2.1\ con- 
verges to the unique solution u to \4-2)^ exponentially as t goes to infinity. More 
precisely, we have the following estimate 

||M(t)-M|p < ||Mo-M|pe-2("+'')* , forallt>0. (4.3) 
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proof. We let v = u{t) in (|4.2p . thus 

i/(Vit, V{u{t) -u)) + a {u, u{t) -u)+ 6(|M|"it, u{t) - u) 
+J{u{t)) - J{u) > (/, u{t) - u) . 

Next, for V = u in (|2.5p . one has 

{u'{t),u - u{t)) + v{Vu{t), V(m - u{t)) ) + a {u{t), u - u{t)) 
+b{\u{t)\°'u{t), u - u{t)) + J{u) - J{u{t)) >{f,u- u{t)) . 

Now, putting together the two previous inequahties yields; 

-{w'{t),w{t))-iy\\w{t)\\^ -a\\w{t)\\^ -b{\u\"u-\u\"u,u-u) >0 , (4.4) 

where w{t) = u{t) — u. From the monotonicity of T{C,) ~ |CI"Ci (14.41) imply 
that 

|||«;(t)|P+2(z. + a)||^(0lP<0, 

from which the announced estimate is readily obtained via Gronwall's lemma. 

□ 
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